The mildly nonlinear imprint of structure on the CMB 
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Abstract 



I outline some nonperturbative relativistic effects that arise from gravitational corrections 
to the Boltzmann equations. These may be important for the study of CMB temperature 
anisotropies, particularly their interpretation. These terms are not included in the canonical 
treatment as they arise from the exact equations. Here a weakly nonlinear investigation of 
these effects is defined and investigated with an emphasis on a Rees-Sciama sourced effect - the 
imprint of structure evolution on the CMB. It is shown that gravitational nonlinearity in the 
weakly nonlinear extension of almost-FLRW temperature anisotropies leads to cancellation on 
small-scales when threading in the Newtonian frame. In the general frame this cancellation does 
not occur. In the context of a flat almost-FLRW CDM model we provide a heuristic argument 



for a nonperturbative small scale correction, due to the Rees-Sciama effect, of not more than 
AT/T ~ 1CT 6 - 1(T 5 near £ ~ 100 - 300. The effect of mild gravitational nonlinearity could be 
more sophisticated than previously expected. 

> : 

: 1 Introduction 

H . 

The 1+3 Lagrangian threading formalism gives generic equations for general CMB temperature 
anisotropies; from these it has been found, using qualitative arguments, that nonlinear terms dom- 
inate on small-scales [22]. In a previous article || I discussed some of the implications with specific 



emphasis on the scattering correction; here I deal exclusively with the scalar gravitational cor- 
rections. The key issue that is dealt with here is that such corrections are not included in the 
canonical treatment^] as they arise from the exact equations; specifically the exact multipole diver- 
gence equations p2[ . The key advantage of using the 1+3 Lagrangian approach for these types of 
calculations is that it is not plagued by the necessity of constructing higher order gauge invariant 
variables when finding nonlinear corrections to the perturbative part of the theory. 
This paper is constructed with the following framework in mind. 

Next, in the introduction, I discuss why such effects could be important. I then summarize 
the general attitude towards nonlinear effects. This is followed by some general comments. The 
main text then follows with section 0, which gives a brief introduction to the 1+3 covariant and 
gauge invariant formulation as reduced to weakly nonlinear form. Here I give an explicit equation 
defining the weakly nonlinear theory in manifestly covariant and gauge invariant form. 

x The linear FRW treatment using the Bardeen perturbation theory 



1 



The core of the paper is a set of calculations using the canonical almost-FLRW results as ansatz 
seed solutions0,so, a brief exposition of the canonical almost-FLRW temperature anisotropy results, 
found using |7|, ||, is given in section |3| This is followed, in Section 3.1, by a recovery of the standard 
derivation of the Rees-Sciama effect; but here in the 1+3 covariant and gauge invariant formalism. 
This will be needed in the last section where we will compare the nonlinear (in the Boltzmann 
equations) corrections arising from the Rees-Sciama effect, to the Rees-Sciama correction itself. 
With the almost-FLRW theory in hand we will then be ready for main section, Section ^. 

In the main section I calculate the gravitational corrections induced via the weakly nonlinear 
corrections using the almost-FLRW theories formulation of the Rees-Sciama (RS) effect over an 
Einstein-de Sitter background (EdS). Section [| is divided into four subsections. The first introduces 
the gravitational corrections in terms of the 1+3 mode functions defined in (Tj, and explicilty 
constructs the Fourier coefficients for the situation of aligned wavenumbers - this simplifies the 
mode-mode coupling considerably. In this section the correction is written in terms of a conformal 
time integration from last scattering until now. It is constructed in a manifestly covariant and 
gauge invariant form for scalar perturbations. In section 4.1 , I calculate the effect in the Newtonian 
frame (see [22, M)f\ - the effect is shown to be suppressed. The third subsection, section 4.2, has 



the calculation in the total frame (see [22, g). The effect is shown to be non-trivial in the total 
frame. The fourth subsection, section 4.3, then attempts to approximate the resulting angular 
correlation function for the total frame effect. This is done by reducing it to a form that is similar 
to the Rees-Sciama calculation given in section 3.1. This then allows me to compare the nonlinear 
correction to the Rees-Sciama one. It all ends with a sparse conclusion and the references. 
As promised, the back-of-the-envelope argument is now given. 

From an order of magnitude argument using the COBE-Copernican limits [ |32| in the exact 
anisotropy equations [^] or Eq. (||) one can approximate the weakly nonlinear correction by 
finding the order of magnitude of the induced temperature anisotropy multipole, ta{- 0(rA e ) ~ 
£O{{a a0 /®)IlA l )- Here we have ignored the acceleration by assuming CDM domination, and have 
used the usual assumption that [f^J < 0|r J 4 £ |; these are consistent with the almost-EGS approach 
[31]. Three results follow. 

First, using ta ( ~ 1CT 5 , for almost-FLRW sources (using {crab/ 9) ~ 10~ 4 ) H and that the 
lower bound on IL^ is that on TA t we find that almost-FLRW sourced nonlinear corrections would 
only dominate the anisotropies at £ ~ 0(|<7 a 6/O| _1 ) ~ 10 , long since diffusion damping and system 
noise have dominated the signal. 

Second, if there are nonlinear sources of shear (due to structure evolution or formation), even 
at only 10 -3 , one can expect an effect near £ ~ 1000. 

An important point here is then that we have entirely ignored feedback of the anisotropy source 
into the correction terms, there could be an accumulation effect. 

Third, we attempt to take the feedback into account by modeling the situation as N repeated 
sources along the line of sight, an accumulated effect : 0(HA e ) ~ VNO(ta £ )- One then finds that 
£ ~ (1/VN) x 10 4 . 

Although this line of argument leads to an unreasonably low value of £, even for a moderate 
number of iterations (or sources along line of sight), one gets a feeling for the situation - the 
universe needs to be very special for there to be no nonlinear effect. We have excluded the effect 
of dissipation via mode-mode coupling, which would smooth the effect and could also push £ up to 
higher values, as will any inherent cancellation effect - as found in the Newtonian frame. However, 
in general, it seems that one should still expect an observable effect. It is for this reason that I 



2 As in previous papers almost-FLRW refers to 1+3 covariant and gauge invariant linearizations while linear-FRW 
refers to the 3+1 linear perturbation theory of Bardeen J5J. 

3 The 1+3 covariant and gauge invariant equivalent of the Newtonian Gauge 
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investigate the gravitational effects in more depth in this paper. 

Perhaps the point is best seen as two fold: first, nonlinear sources of shear (or &h) would 
result in a moderately amplified correction apparently not accounted for in the canonical treatment 
because of its, by construction, linearity, second, nonlinear feedback (or, if that does not do the 
trick, accumulated sources) could have an even more noticeable effect. 

The general attitude towards the importance of nonlinearity in both observational cosmology 
and CMB studies can be understood in the context of popular lore surrounding the Rees-Sciama 



effect (see for example [23, 

The two points that are usually made with regards to the Rees-Sciama effect are: first, that 
numerical simulations (where nonlinear effects dominate the linear effects only above £ > 5000 or 
that the effects are a problem near to a sensitivity of 10 -7 ) corroborate the scaling arguments and 
second, that the scaling arguments indicate that the effects are ignorable: that the effect suffers 
cancellation as (kdrj)^ 1 for 5r] being the time-scale of change in the potential |fT6| l. The idea is 
that the Poisson equation relates potentials to densities as (krj) -4 while the volume from the mode 
coupling introduces terms like k 3 meaning that the effect scales as k 2 P{k) - which is ignorable, 
particularly when compared with the dominant Vishniac corrections (which cannot be ignored). 

The novelity of using the Rees-Sciama effect as the source of nonlinear gravitational effects is 
that the peak in the nonlinear correction seems to be near to the peak in the radiation transfer 
function, I ~ k(j]Q — ?/*), somewhere near £ = 100 — 300 [3C[. The additional £ scaling from the 
nonlinear correction then introduces an additional k 2 scaling in the angular correlation function. 
It is for this reason that the scalar gravitational correction will be investigated in substantially 
more depth using the Rees-Sciama corrections as the source term. It does not seem unreasonable 
to suspect that it could dominate the acoustic effect. 

The additional linear scaling in I ( p^l and Eq. (||)) may cause the scaling argument to fail, 
nonlinear effects could become important well before £ > 5000. It is the latter point that needs to 
be properly understood, especially given that this scaling does not seem to arise in the perturbative 
analysis; it appears as a nonperturbative feature of the exact small scale treatment. 

In addition, many established treatments carry out the analysis in the linear gauges which set 
the shear to zero : such as the longitudinal, Newtonian or conformal Newtonian choices [27, 23 1, 
which have no nonlinear corrections of the form discussed here - this will be shown here in detail. 

A second order EdS extension of Pine-Carrol [^7|] was carried out by Mollerach and Matarrese 
[26] using the Poisson gauge, which gives the appropriate Newtonian extension at higher order; 
it requires that D a <7 a 6 0, somewhat less restrictive than a a b « 0, it is however a perturbative 
analysis and so will not uncover the corrections discussed here. The main problem with a strict 
perturbation theory approach is the need to maintain gauge invariance at each order [^, gq]. In 
addition the strict perturbative approach is a computational tour de force, even in the case of 
investigating "mildly nonlinear effects". 

The points that are being made here are: (1) that when calculating the influence of nonlinearity 
on temperature anisotropics using a Newtonian gauge choice within the linear theory context one 
should view the results of such calculations with suspicion. This is because: (i) the gauge is 
inconsistent at higher order [33|, and (ii), suppresses the mode- mode coupling artificially - the 
cancellation should not be a suprise. (2) When using a higher order perturbation theory, one 
should not expect to find these effects either; the effects discussed here are not of a perturbative 
nature. 

Although these arguments may seem esoteric and somewhat obscure, the point at hand is 
understanding the stability and validity of the use of the almost-FLRW models on observation- 
ally interesting scales. Can the observations be explained using non-almost-FLRW corrections? 
are non-almost-FLRW corrections something worth worrying about? From the perspective put 
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forward here: the resolution of this is not obvious given the nature of the construction of the 
canonical linear-FRW programme. The situation may become even more complex if one has to 
worry about a gravitational wave (GW) background (or other large scale effects). Perhaps another 
way of phrasing the question is: why should the time derivatives and spatial gradients of the dy- 
namical and kinematical quantities be small?, when realistic cosmological models are apparently 
stable to inhomogeneity, particularly given the subtle nature of realistic models which have coupled 
inhomogeneity and anisotropy. 

I now define the weakly nonlinear theory. 



2 Weakly nonlinear almost-FLRW 

The almost-FLRW model is expressed in terms of the 1+3 covariant and gauge invariant pertur- 
bation theory. The temperature anisotropies are defined in terms of the covariant total direc- 
tion brightness temperature and the bolometric average in the almost-FLRW universe: T(x, e) = 
T(x)(l + t(x, e)) |22|, ||]. The almost-FLRW temperature anisotropies for generic-^ (in particular 
£ > 2) were the focus of § while the implication of the observational constraints on the dipole, 
quadrupole and octupole were the focus of 1^] - where the theoretical foundations of the use of the 
postdecoupling CDM dominated almost-FLRW models were established, and hence the plausibility 
of the use of almost-FLRW models on large scales within the context of relativistic cosmology. 

The notation used here follows [|], as in the previous papers; in summary the temperature 
anisotropy r(x,e) has been expanded in terms of a multipole expansion, r = J2e T A e O Ae with 
Ae = e^ A ^ forming a Projected Symmetric Trace- Free (PSTF) basis 0]. The equations are all 
expressed in terms of the 1+3 threading variables derived with respect to a preferred u a -frame @, |5|: 
Vfe« a = ^habQ+o-ab+eabcUJ — A a Ub, defining the expansion: D a u a = 0, the shear: "Duu^s = a a b, the 
vorticity vector: oj a = — ^e a bcD b u c and the acceleration: A a = u a = u b VbU a . The dot is indicative 
of comoving time with respect to the u a -frame, while a prime, will be used to denote conformal 
time. The spatially projected covariant derivative is denoted: h ab Vb = D a . 

Here u a u a = —1, and the direction vectors e a are defined such that e a e a = +1 and e a u a = 0: 
the momentum is the given by p a = Eu a + Xe a such the E 2 = m? + A 2 . We refer the reader to 
references Jl2| , |], ^2|, [B| for further details on the formalism - though care should be taken between 
the differences in the formalism in |J and [12, ||, 22], I will be following the later. 



The Weakly nonlinear equations for the brightness temperature anisotropies are defined by: 
f + e a D a r + C[t] + (5f) ov + (5r) NL + SC N 

almost-FLRW weakly nonlinear terms 

Here we have used the exact formulation found in |2^] in the small scale (high-£) limit and FLRW 
linearization discussed in to construct the weakly nonlinear theory. The respective terms used 
in equation (Q) are, the linear order coupling between the anisotropies and the field equations, 
the linear order Thompson scattering correction, the Ostriker-Vishniac (OV) correction and the 
dominant nonlinear (NL) corrections: 

B{x t e) « -lB a T a + (D a lnT + A a )e a + a ab e a e b , (2) 

C[r](x,e) « -a T n e (fie a v* + f 2 e a e b r ab - t), (3) 

6t ov ps -a T n e p ab B a pv b , (4) 



1 ov 
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_( A (a tT A e _ 1 ) _ l A{jT bA^ + u b ebc {a tT A e _ lC )l j (5) 

5C NL « +0 - T n e O A Jr^-i^ ) + ir^X B l • ( 6 ) 



As was shown in [22|, the dominant nonlinear contribution to the temperature anisotropics at 
high-£ will arise from term (Eq. ||) above; if gravitational nonlinearity between decoupling and now 
significantly alters the small scale angular correlations the dominant contribution will arise from 
this term. By considering this term in isolation, that is, outside of the context of the exact equations 
and within the context of corrections to the almost-FLRW theory I am trying to put bounds on 
the possible effects that this term may have on the small scale angular correlation functions and 
hence on the acoustic peak. 

One notices that for almost constant coupling one may expect to find an effect that increases 
linearly in I. Furthermore, these equations are valid for the so-called scalar, vector and tensor 
perturbation split - at least near to linear order where such a split is meaningful in the almost- 
FLRW theory. 

The scattering term includes two weighting factors which are used to deal with the anisotropic 
scattering correction and the polarization correction; this has been discussed in detail elsewhere 

0,0, ©a- 

The third and fifth terms are not dealt with any further in this paper, that is the Ostriker- 
Vishniac (OV) correction and an additional nonlinear scattering correction. The OV contribution 
represents additional variations in the baryon streaming velocities, induced by gradients in the 
matter energy flux that generate flows perpendicular to the line of sight e a . 

A last comment on other outstanding effects is with regards to the Sunyaev-Zel'dovich (SZ) 
effect, in particular the kinetic-SZ effect, this arises through e a v^ in re-ionized re-scattering. The 
thermal-SZ effect arises from including a electron-baryon pressure term. The kinetic effect is 
accumulative, and is due to along- line-of-sight streaming velocity differences - this tends to cancel 
in much the same manner as the baryon infall effect does during slow decoupling. 

In order to maintain our focus here: we are primarily interested in bettering the understanding 
of the gravitational nonlinear effects on the CMB temperature anisotropics (Eq. ||). This is to 
attain two immediate goals: first, to demonstrate how these effects would arise as a correction to 
the standard gauge invariant treatment following the canonical treatment, the work of Hu et al 



[13, 14] in the formalism of Bardeen [jlj] and Wilson ]p5| , |34j], and second, to demonstrate how to 
implicitly calculate such corrections analytically by finding the corrections to the mode coefficients 
of the temperature anisotropies. 

I will first investigate the scalar effect in the conformal Newtonian frame using weak-coupling 
[15] and high-^ approximation techniques (see latter sections and p2|). These calculations represent 



new effects that do not arise in the canonical treatment. As pointed out before, a caveat here is 
that in the Newtonian frame any extension to second order must be treated with extreme care. 
Such a threading is inconsistent beyond linear order H - we are relying on weak coupling and 
sufficiently small peculiar velocity corrections in order to argue validity and hence consistency of 
our treatment. 



4 The energy flux in the Newtonian frame with respect to the total frame is : Db<ja ~ D(,g a — (D^p + DbP)v a - In 
the situation with dominate pressure free dust one then has an additional scattering contribution as : 



<5T ov (r?o) « - / («'e K ){D a p B v a }dr). 



Here, the divergence of the energy flux will vanish along the line of sight by cancellation due to falling in and out of 
potential wells on small scales - however the tranverse components of the quantities in the brackets will not vanish 
generically. 
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Ideally such a calculation should be carried out in a manifestly gauge invariant manner (using 
the generic Lagrangian threading in the total frame) or in frames that are known to be consistent 
at higher order; such as the constant expansion frame or CDM frame §: we do this for a nonlinear 
Rees-Sciama imprint on the temperature anisotropics. With such frame choices one then has the 
additional freedom of moving away from weak nonlinearity. We now begin the explicit calculation 
of the weakly nonlinear correction (with no nonlinear feedback beyond the effect itself) using the 
weak coupling approximation [jD] ]. 



3 The almost-FLRW anisotropy sources 



The calculation in this section aims to, first, provide the background so as to, second, reproduce 
the results of Martinez-Gonzalez et al [|4|] and Seljak [pcj]. 

I use the temperature anisotropy results for almost-FLRW models as found in ||. The key 
worry in this calculation is the possibility of cancellation between the primary source feedback 
and the secondary source feedback. The almost-FLRW primary (P) source || |l3[ (using matter 
domination to get ~ — t/q is the conformal time now and rj* is that near last scattering) is: 



{21 + 1) 



[5T-<P H ]{k,rj*)j e {k{ m -rj,)) 



The almost-FLRW secondary (S) source [y, 13 1 is: 

Using weak-coupling (W) and not specifying the end time (7/0 — > rj) the latter becomes 

A21 + 1) rru, 



t-; w (m) 



Pi 



.^^> H (k, m )e- K , and m « 77 - (2£ + l)/2k. 



(7) 



(8) 



0) 



To summarize the power spectrum definitions^, recall that: P(k) = | A(/c, ?7o)| 2 , for A = D(rj)A(k, 770) 
and D(rf) ~ (t/ 2 /??o)) m addition^ one has 



|ff 2 O A(?7,fc) « -k 2 (a<S> A (k,0)) 
I can then write the primary and secondary sources as (cf. |l 



Tj{k,r]) 



r! w (k,v) 



(21 + 1) 

Pi 

{21 + 1 

Pi 



5T --Hfa k- 2 (-^)A(k,7 l0 



D t 



j e (k(rj - 77*)), 



* 1 „_« ^ 



A(fc,7? ) 



(10) 

(11) 
(12) 



I assume the standard Einstein-de-Sitter (EdS) results: D (7/) = a(t)/ao from a/ao = t? 2 and that 
P(fc,t) = P(k,0)D(t) along with the useful relation {Q D*/a*) 2 ~ ft 1 ' 54 @. The intention here 

5 Note that the definition using 8 k differs from that using A by fc 3 : 

A 2 = ^cxfc 3 P(fc) for, V(k) = (\5 k 5 k> \ 2 ). 
a In k 

I will be using the dimensionless form; the variance per In k. 

6 ±D a lnp B w p-^Eab D a l n/ 9 B « +f A(fc,?7)Q a and £„6 « iD (a D b) - ^(x)) « +D (a D i , ) $ A (a;) 

=> |i7 2 fioA(fc,77) 



-fc 2 $ A (fc,0). 
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is to use the EdS growth factors with a nonlinear correction to A (A;). As pointed out before, we 
are modeling the effect of small-scale nonlinearity in the matter density so as to recover the Rees- 
Sciama like corrections for the exact theory. In addition one could include damping as k! /k to 
include cancellation and diffusion (see jO]]). 

It seems reasonable to assume that the dominate contribution will arise from the ISW terms 
via a Rees-Sciama correction; this will dominate the early-ISW, Sachs- Wolfe and acoustic sourced 
primary effects; assuming that there was little or no nonlinearity near last scattering. So I will 
only consider the secondary sourced effect. 

3.1 The Rees-Sciama effect (RS) 

In order to understand the nature of the nonperturbative corrections we need to understand the 
Rees-Sciama effect. Here, as promised, I reproduce the well know result using the results from 
the previous section. From the almost-FLRW temperature anisotropies, using that <£#(&;, 77)' ~ 
D'(k,r])<fr(k,r]o) [] and following [ 30 1 , but using the covariant and gauge invariant notation (Eqs. |9| 
and|): 

rf(k,ri ) « -2^^^^e- K D'(k, m )^H(k,Vo), for k m >£. (13) 
This term is zero for krjQ < I. Now using the definition of the angular correlation function M: 

Pt / t.2 j Tj ^RS/u m \|2 



c r = ^^[fwL kdkVi " (Mo)| ' (14) 

We find the angular correlation function and the Rees-Sciama effect in terms of the power spectrum 
of the potential: 

4 r°° 1 
Cf «-(4vr) 2 y o k 2 dk(D>(k, m )) 2 Pz H (k)^. (15) 

The power spectrum of the time changing potential is defined as 

Pv H (k, V ) = (D'(k,r ] )) 2 P* H (k) (16) 
[30] to find that (here we are using r\i = rj — tjk such that k = £/(t]q — rj) for all kr) > £): 

A roo 

Rees-Sciama: Cf « -{4tt) 2 / dkP&Jk^o - f). (17) 

I can rescale the integration variable from k to rj, dk ~ \£/{f]o — rj) 2 \(dr]) : and change the integration 
limits (the integral is non- vanishing for k > i/rjo) to write the angular correlation function in terms 
of the area distance r « (770 — i]) : 

«4(4vr) 2 r ^Pv H (£/r, V ). (18) 

J * 

This recovers the result of |M]. From |3^] we have that 

P *>n = \% a ') 2p ^ k )- ( 19 ) 



7 This is nonlinear for the CDM dominate flat almost-FLRW models - where $H(k,vY ~ ®ii(k, 0)' w 0. To 
understand the notation used here. The Rees-Sciama correction arises from (i) the term < J > £r(fc, r/) — D(k, r])&H(k, 770), 
while the (ii) generic ISW effect arises from <E>£r(fc,0)' which is written as (D(rj)^H(k, 0)/a(rj))'. 
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This then gives the standard Rees-Sciama effect. The nonlinear power spectrum is approximated 



from [24, pQj : P$i (k) oc k 4 -P( 2 )(/c). The relationship between the second order power spectrum 



and the linear one can be found from 

roo 



roo 

P(2)(k)^J o q 2 dqP(k)P(\k-q\)D 2 {2) (k,\k-q\) (20) 

where D(2) is the second order growth factor, and P(k) = A{k/ko + k 4 )~ 1 e~ ka for ko = /i/30Mpc, 
a = 8/i -1 Mpc such that for k < ko one has P(k) oc k and for a -1 > k > ko one finds P(k) oc 1/k. 
In this latter region one is lead to find: P(2)(k) oc k^P(k). With the additional cancellation as 
(k5rj)~ l one finds the usual result that the source scales as k 2 P(k) in this region of interest - as 
pointed out in the introduction. Here one then finds P^ A (k) cx k~ 2 P(k) as the power spectrum of 
the time rate of change of the potential in the region a -1 > k > ko, leading to the usual scaling 
argument. 

We are now ready to investigate the new effects in more detail. 



4 Weakly nonlinear gravitational corrections 

I use the usual scalar mode decomposition but now include a mode-mode coupling between the 
wavenumber (assuming that the directions are all aligned along the line of sight) 8 . 

The nonlinear corrections multipole coefficient and the temperature anisotropy multipole coef- 
ficients can be put into mode coefficient form: 

Kl = E(^)^< and =H^Qt- (2i) 

k a k a 

Using the Fourier convolution theorem, with a little algebra and integrating over the solid angle, 
these, along with the mode form of the scalar potential ||, using the PSTF relation e ae+1 Ae+1 
0] and an identity for the delta function, 5(2k) = ^S(k), will allow us to write Eq. (||) as a mode 
coefficient for the correction in terms of those for the temperature anisotropy and scalar potential. 

Now there are four different cases; using in addition that I 3> 1 and that e a O aAe = (£ + l)/(2£ + 
1)0**: 

roo 1,2^1, 

A { a e r Ae _ l} ~ (4n)Q Ai (x a ,k°) J q —A^n.^k* - fc|), (22) 

i 1 roo U^fjU 

A a r aAe ~ ^)-^—Q Ai { x -X) J o j^A(k)r £+1 (\k* - k\), (23) 
" { a,a,. l rA,_,)~(^)QA,(' a ,K)l p^(fc)T,_ 2 (|k* - k|), (24) 

°*>r Ml ~ (4,)0^(xr fe ;) ^±M±| 1°° ^(k)r M (W - k\). (25) 
This approximation scheme only holds for high-£, this will then give 

roo b^rlU 

A {ae TA^) ~ {^)Q Al {x a X) _A(fc)T*_i(|A:*-fc|), (26) 

8 We do not need to worry about the direction couplings then (which is of course essential for the Vishniac effect 
and the nonlinear scattering correction which would generate a nonlinear Vishniac effect). We therefore do not have 
terms such as: \ £ fe „ (A(k a )M(k* a - k a ) + A{k* a - k a )M{k a )). 
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A a T a A f 



1 roo hZflU 

(An)-Q Al (x a ,K) —A(k)r e+1 (\k*-k\), 



a 



(aeae-i T Ae_ 2 ) 



(^)Q Ae (x a ,K) 



<y ah r ahAl ~ (4ir)Q Ae (x a ,K)\J k 



k l dk 
■° k 2 dk 



k (2vr)3 



*{k)T t - 2 (\k* ~ k\), 

a(k)r e+2 (\k* -k\). 



(27) 
(28) 
(29) 



I drop the vorticity (we are assuming matter domination for which it becomes natural to treat 
vorticity effects as very small). I then have the weakly nonlinear correction on small scales, Eq. 
(p(), in mode form: 



-(4ir)£Q Ae (x a ,K) 



k 2 dk 
k W? 



a(k,t)T £ _ 2 (\k* - k\,t) + — a (k,t)r e+2 (\k* - k\,t) 



- A(M)7*_i(|fc* " k\,t) + -A(k,t)T £+1 (\k* - k\,t) 



(30) 



I am interested in finding the effect on scales : k > k$. The above equations are valid for any 
it a -frame as I have not yet frame-fixed the theory. On dropping the mode functions I find the 
following equation: 

roo ip'AU r 1 

' * ' a(k,t)r^ 2 (\k* - k\,t) + — a(k,t)ri +2 (\k* - k\, t) 



d 



dt 



(6r)^(k*,t) 



k (2tt) 3 
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1 



- A(M)T*_i(|fc* - k\,t) + -A{k,t)r e+1 (\k* - k\,t) 



(31) 



Finally, we can write this out in terms of the conformal time through: dt = adr/, and write the 
effect today in terms of a timelike integration (under the assumption that the effects are secondary 
- as opposed to primary): 

pr] poo hZrfh. 

5T^{k*, m )K-(A-K)ij^ aWdnj —-[a(k lV )r e - 2 (\k*-k\,ri) 



+ ^-a(k,t)T £+2 (\k* - k\,rj) - A{k,rj)T^i(\k* - k\,rf) + -A{k,rf)T t+1 {\k* - k\,r/) 
lb 4 



(32) 



The question is: "How does this compare to the linear primary and secondary sources from the 
canonical treatment? 1 

Towards answering this two distinct effects need to be highlighted: the first, is the result of 
mode-mode coupling between the linear primary and linear secondary anisotropies (Eq. |l2|) and 
the linear shear and linear acceleration sources; this will most likely lead to a small smoothing 
effect at high-^. Second, the effect of local nonlinear matter dynamics coupling through into the 
anisotropies. It is this latter effect that I will emphasize here, because it promises to be dominant. 

4.1 The Newtonian-frame correction 

Once again starting from Eq. (||), the high-£ correction assumes the form 

(6t) a * « -£ {jO hc T hcA t + o-<«/«i-ir^-2> - A^t a ^ + \A h r hA *- + cA 6f >r^- lC >) . (33) 
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in the Newtonian frameQ, where equation (j33|) can be written in terms of the acceleration alone @. 
Then it can be written, for scalar perturbations, in terms of the Newtonian potential f|: 



1 7 



bA, 



Here I have used that A(k) 



i 1 (jD< a $ A )r A 
k a \/a)$ A (k,r]) for A a « D a $ A . 



1 



(34) 



4.1.1 Expanding out the correction 

I change from comoving time, t, to conformal time, 77, in Eq. (|34|). Invert the resulting equation 
into an integral equation and now work in the conformal Newtonian frame using the conformal 
time variable, dt = adr/. The integral inversion will only work for weak nonlinearity as we have 
excluded the general feedback: 



(Sf)^{k*, V ) « -{4tt)1 



* r 00 k' 2 dk' k! 



'v* Jk (27r) 2 a 

It becomes convenient to rewrite this in terms of the following sources: 



n-i{k - k',rj) - jTt +1 (k - k',r)) 



T e (k*,r,)=Tl(k*,r,) + Ti(k*,r l ) 



T e -i(k*,r)) - -T i+1 (k*,i]) 



(35) 



(36) 



t v " (2£ + l 

where S and P denote the secondary and primary induced corrections. 

Essentially the argument in the following two sections is that the time integration over Tn{k,r]) 
cancels when using weak coupling on small scales, that is for high-£ 



4.1.2 Primary Sourced correction 

Using the almost-FLRW integral solution || for the temperature anisotropies in the mode coefficient 
formulation of the nonlinear correction (Eq. |35j) and using (|36[) : 

^§^(V0,k*)^+2£ J"°d V J™^k<t> A (k, V )Tf(k*, V ), (37) 



(2£ - - fcj) _ 1 (2£ + 3)j e+1 (\k* 

4 



(38) 



4 f3 i+1 

Notationally, what is important to realize is that the Bessel functions have conformal time argu- 
ments as tjq — r] inside the time integration. The co-efficients inside the square brackets are evaluated 
at 77* while the ones outside at r/. The best way to think of this is that the integral solution is 
integrated up to some arbitrary time parameter 77, this is then convolved as in the integral above 
and integrated from 77* until now, 770- 

Now, I will again use the weak-coupling approximation [15|, arguing that the variations in the 
CMB temperature anisotropies are significantly more rapid than those in the potentials. This 
will allow me to reduce the potential term and remove the conformal time integration down the 
worldline: 

no fW k 



(21 + 1) 



+ 



2 \k* - k\ 



[6T + $ A ](\k*-k\,ri.) 



X 



(21 - l)P t 



(39) 



9 From Iql, where u a — n a w u a + v a such that D( a 7i(,) = a a b = 0, which if dominated by pressure free dust leads 
one to expect uj a ~ 
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We should expect some cancellations to occur between the k* and —k* terms too. I use that 
= £/{2£ - 1) and Pe/Pe+i = (2£ + !)/{£ + 1) 0. This allows the following reduction: 



(21 + 3)0 e 



+ 3){2£ + l) 

3 
I 2 



and 



(2* - l)ft 



I now take the high-£ approximation (large- A;), again: 

1 + 3)/(£+ 1)1 w2V£ and ^/(2^ + l)\A^" 1 ) 



(40) 



(41) 



to find that by putting Eq. (M) into Eq. d39| ) and using Eq. (|4l|), that upon collecting terms to 
cancel factors of 1/4, we recover : 



0i 



{21 + 1) 



l In 



d V k$ A (k, V )T[(k*,ri) 



i nr 



k\ 



(42) 



:[6T + $ A ] ($ A (k, v f_- k ) - $ A (k, 



v k Zx k : 



Putting this back together with Eqs. (37 and |42| ) I find (the additional £ in front of (|37| ) leads to 
the \f£ scaling below and the factor 2 cancels): 



(2£ + 1) 



ivo,k* 



00 k 2 dk k 
k (2vr) 2 \k*-k\ 

:[5T + ^ A ] (<S> A (k, V fSx k ) 



k*-k- 



(43) 



First, for sufficiently high-£ (that is small scales) we could argue that the peaks of the Bessel 
functions are sufficiently close together (again) such that % ~ %±i; one immediately notices that 
all the relevant terms will cancel out in Eq. fl43| ) , we then find that the primary sourced corrections 
vanishes 



(2£ + l 



-(«5r)r- p (%,fc)~0. 



(44) 



So it would appear that in the weakly nonlinear case there is no effect. The cancellation has nothing 
to do with the damping factor which we can safely ignore. 

I have shown that the primary sourced effect cancels (this I have labeled NL-P). I now show that 
the secondary sourced effect cancels too. This is more important as this is how the Rees-Sciama 
effect would arise in the Newtonian frame. 



4.1.3 Secondary source correction 



01 



From Eq. (35) the term of interest is 

Tt(k',r,)= {n+l 
Now from the ISW source of anisotropy I have Eqs. (|l2| and 



-(21 + 1) 

01 



k* 3 



sw(% ! ' 



k*- 



3Hn 



D /£>' 



0"0 



a \D 



(45) 

(46) 
(47) 
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Putting Eq. @ into Eq. @ I find on using that fie/Pe-i ~ \ and that Pe/Pe+i 

21 p t \ nrfA(k*, Vo y 



2 (for 1): 



+ 



./%_! 2^; V 2^ V fc* 3 
1 / 2^ /3A fA(k*,rj ) 



4 VA+i 2^7 V 2£ V A;* 3 



for 

ISW (Ci) 



'iswV'«-l 

SI 



lA(k*,r] ) pR fa ; i ,* 

a iswl%-l 



'^(S ISW (rfc 



Now, using that « % +1 



2 A;* 3 V 2^ 
iji ~ ij — £/k, I then find that 



(48) 
(49) 

(50) 



I have arrived at advertised result; that acceleration sourced nonlinearity is at best small and in 
the limit vanishes. 

This means that in the Newtonian frame there may be mild nonlinearity (in the sense of second 
order perturbation theory as pertaining to the dynamics of the matter) but there are no weakly 
nonlinear effects due to the acceleration potential in the Newtonian frame - as arising from non- 
perturbative small-scale effects in the radiation dynamics, which are the dominant gravitational 



effects on small scales [22 



The important point is that this provides some indications that the effects of weak nonlinearity 
are suppressed in the Newtonian frame; in some sense the Newtonian frame calculation is stable 
to weakly nonlinear contamination making it an ideal first approximation when trying to include 
additional nonperturbative source terms. 

In the Newtonian frame one needs only include the linear Rees-Sciama effect (Eq. [17]) which is 
known to be small |30f l. In the CDM- frame (or the total frame) one would need to be more careful 
with the nonlinear Rees-Sciama corrections (NLRS). 

4.2 Total-frame correction 

In the CDM dominated almost-EdS universe we have that A a ~ 0. One then only needs to consider 
the nonlinear corrections due to the shear (given that the vorticity is small) - this is a nonlinear 
shear sourced correction. 

Here I consider the nonlinear corrections (Eq. |33| ) due to the shear in the generic u a -frame (Eq. 
32 in its manifestly covariant and gauge invariant form): 



(5Tf £ L (k\rio)^^£ dr,' 



no 



k 2 dk 
k W? 



a(k, rj) 



n„ 2 {\k*-k\) + ±-T l+2 {\k*-k\) 

lb 



(51) 



This will not generically vanish on small-scales; however in the case of sufficient matter domination 
I would expect the effect to be small. In the context of an EdS background the ISW sourced effect 
will very small, because the effect of {<&H{k,rj) / a)' « (D/a)(D' / D — a'/a)&H (k, 0) is vanishingly 
small - which sources the ISW terms in the temperature anisotropy. So I need to consider the 
nonlinear-ISW effect - the Rees-Sciama effect. 

The coupling between the primary sourced effect and the linear-FLRW shear could be expected 
to, at best, lead to a smoothing of the acoustic peaks; this still needs to be properly shown. If we 
interpret the shear as sourced by nonlinear dynamics (which leads to the Rees-Sciama effect - the 
imprint of nonlinear dynamics on the temperature anisotropics) I can then isolate three possible 
effects: 
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1. The primary correction about the flat CDM background: due to the coupling between the 
almost-FLRW shear (a term like (a' /a)$A(k, 0)) to the primary sources of anisotropies - the 
acoustic and Sachs- Wolfe effects. This may be important on intermediate scales. 

2. The primary nonlinear correction about a flat CDM background: due to the coupling 
between the primary sources of CMB temperature anistropies and nonlinear shear, induced 
by nonlinear small scale matter dynamics. This would be important on small scales. 

3. The secondary nonlinear correction (or nonlinear Rees-Sciama effect) about the flat CDM 
background: due to the coupling between the nonlinear-ISW effect (Rees-Sciama effect) in 
the CMB temperature anisotropies and the nonlinear corrections to the shear (a term like 
(D' /a)§A(k, 0)). This effect is expected to be the dominant small scale correction - and 
is due to small scale nonlinearity in the matter dynamics generating radiation anisotropies 
which couple back to the nonlinear shear. 

I focus on this latter possibility, with the intent of getting a feel for the form of the angular 
correlation function due to this effect. The interest in the secondary correction, a Nonlinear Rees- 
Sciama effect, lies in the realization that the peak in the Rees-Sciama effect is in the range of 
£ = 200 — 300; where one expects to find the peak in the angular correlation function (apparently 
due to the projection of Doppler and acoustic oscillations pSf] ). I now provide a simple calculation 
demonstrating this contribution. Once again this calculation should be viewed as the departure 
point towards more sophisticated treatments. 



4.2.1 The almost-FLRW shear source 

The shear in the energy-frame (which coincides with the total-frame in the CDM dominated case 
considered here), is given from || as: 

\{p + p)<Jab « (D (a D b) $ H )' + 3HD {a D b) $ H ~ HD {a D b) (<S> H + $ A ). (52) 
In addition, for sufficient matter domination || (A a ~ and hence &h ~ — §a) : 

pivoWb ~ -2(a 3 D (a D 6) $ H )'. (53) 

Now a 2 p(r]o) « 3KQo m /Qo m — 1 ~ 3ci 2 HqQo (where K is the FRW curvature) and p(rj) = a~ 3 p(r]o), 
which for scalar perturbations J?], ^| gives the mode coefficient for the sheaif°| using a flat CDM 
dominated model: 

3 77 2 0n / k 2 \ ( k 2 \ k 2 

-^j^a(k, V )Q ab « -—<S> H Qab )+?>h[ -—§ H Qab ~ --H$ H (k, 0)Q ab . (54) 

I then find, on using that rj) = D(r], k)^n{.k, r/ ) and $h(A;, r/o) = &ii(k, 0), here for the dust 

case, and along with a change to conformal time {dt = adn), to find 

3H 2 n . . k 2 . . 

--^-a (k,rj) « --—$ H (k,r, ). 55 

I have in mind ^> H^h for <&#(&, rj) = D(i])&jj(k, 0); the small scale nonlinear situation, where 
in the linear case D = a. So I am then able to use the following result: 

~H 2 n o-(k, V ) « -k% *h(A;,0). (56) 

10 Using D a $H = +^H(k,rj)Qa and D {a T> a) $ H = -^®H(k,7i)Q ab . 
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Equation (j56|) gives the almost-FLRW shear in terms of the scalar potentials 0) ~ — $#(&, 0). 

However, more importantly, it also shows how the nonlinear dynamics will affect the shear in the 
1+3 covariant and gauge invariant approach, via the growth factor D(k,rj). The nonlinear source 
of <3?# + H^h will give a mildly nonlinear shear effect. 

An important point here is that we have used the shear in its total-frame formulation (that 
natural to the Lagrangian threading), while the CMB temperature anisotropies have been found 
in the Newtonian frame (that given by an Eulerian threading in the exact theory) . To understand 
why we can use these two together we recall (i) that E a ^ is invariant under linear order frame boosts 
||, (ii) Tt (k, rf) is invariant to linear order for I > 1 (we already have the source terms for the CMB 
temperature perturbation, 5T(k,n), and the dipole, Ti(k,n), from || in the Newtonian frame). We 
can use the scalar sourced anisotropies as found in the Newtonian frame || - we use the shear in 
the total frame along with the evolution equations for the other various dynamic quantities^. 

Notice that from the div-E equation, k 2 ^n(k, 0) m +^(D/a)H^Q A(k, rj ) ||. Now I have that 
A(fc, 77) ~ a5 + a?52, where in an EdS spacetime only the second-order part contributes to the time 
changing potential, hence to the growth factor as D ~ D + (j]) ~ a?{rj). So I then find: 



D'(V) 
a(rj) 



Anl(Mo) 



(57) 



which is for sufficient matter domination on small scales. 

In this section I have derived the form of the shear source term, I will use this result in the next 
section: cr NL ~ a. 



4.2.2 Nonlinear Rees-Sciama effect 

The nonlinear Rees-Sciama effect will take the form 



21 + 1) 

gNLRS 



, i * \ ~ & (a o\ r a r k2dk «ni 



a(k, V ) [r^ 2 (\k* - k\,n) + ±TfS 2 (\k* - k\, V ) 

with the linear Rees-Sciama source terms 

D' 



a(k,n) 



Tl s (k,n) 



-A(fc,7? ) 
(2£ + l) 



1 



-k z $ H (k,0) 



D' 



tt e 



21 k 



S '0 

D>( m )<t> H (k,0), Vk v >£, 



(58) 



(59) 



(60) 



where % ps 77 — (£ + \)/k, as before. Now substituting Eq. ( ^9|) and Eq. ( pC|) into Eq. (|58|) I am 
able to find: 



,NLRS 



(21 + 1) 



-{k*,7]o) 



8tt 



o S2 o 

k 2 dk 
k Wf 



-(e" K ) 



\k* - k\ 



(61) 



$ H (k)$ H (\k*-k\) 



U\k*-k\), 



lA\k* 



k\) 



w D'(ri) , 
v* a (V) 



\k*-k\ 



(62) 



n The Einstein field equations are solved in the total frame. See Q for the scalar theory, the tensor part is given 
in |l2|]). The scalar theory is solved in terms of the scalar potentials, §A{k,r]) and ??)), and a peculiar velocity, 

v(k,rf). These can then be substituted into the Newtonian frame solution 
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where I have used the high-^ assumption, which gives the ratios of the coefficients 

A? 



4' 



and, 



01+2 



~ 4, W > 1, 



Using that r^_2 ~ %+2j I can show that : 

- k\) ~ i6 



5 2r3P 



16 



1 _ 1 1 

5 2 |fe*-fc| 



(63) 



(64) 



I also need to use that the second order growth parameter is -D(a) ~ a 2 (??), and that for £ ~ fc, the 
higher order terms will drop off faster under the fc-space integration. 



4.3 Approximating the angular correlation function 

In order to find the angular correlation function I first find the mean-square of the correction to 
the CMB temperature anisotropy. 
I use thatQ 



(* H (k* - k)<S> H {k)<S> H {k" - k')<f> H (k')) = (27rfP t , H (k* - k)P^ H {k) 
5(k* - k") 5(k - k') 5(k* - k") 5(k* - k - k') 



k* 2 



k 2 



+ 



k* 2 



(k* -kf 



(65) 



to then find the mean-square of the correction to the CMB temperature anisotropy due to the 
nonlinear Rees-Sciama coupling between the local matter nonlinearities and the CMB temperature 
anisotropies. 

It is worthwhile pointing out that the basis of the argument used to construct the maximum 
mean square contribution is that it will be near the peak in the radiation transfer function: k ~ 
\k* -k\~ £(rj - rj) This, along with {D') 2 P^ fs P^,, shows 



(21 + 1) 



\5rf hRS (k*,rj))\' 



(4tt) 5 



8tt 



o S2 o 



(66) 



J k k 2 dkp^ H (k, m )p^ H (k* - k, ve )- 



2 \k* — k\ : 



The correction to the angular correlation function due to the nonlinear Rees-Sciama terms Eq. (j36|) 
is constructed from the definition of the angular correlation function [Q, ||]: 



£NLRS 



ft 



k* 2 dk*\5^ s {k*, m )Y 



7T {21+ l) 2 

I then find that the approximate maximum correction to the angular correlation function is: 
8vr 



(67) 



£NLRS 



16vr z 



1 2 



3H2n 



dk* 



fc 



k dkP& H (k, r]o)P$> H {k* - k, %) 



k2 6(k*-£r, ) 



oj 



12 This arises from 8(k a - k' a ) = ^8{k - k')S(e a - e' a ). Recall that (|$(fc, 0)$(fc', 0)|) = (27r) 3 P $ {k)S(k a - k' a ) 
1 To many it may seem as though this is being written in by hand - this is not the case. One expects the Rees- 
Sciama effect to peak between, i = 100 — 300 pOj, while the peak in the radiation transfer function is also in this 
region in the standard CDM model M . 
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I have approximated the time integral by using the assumption that the effect is most prevalent 
near t ~ k(r/o — rj) (again). We use that P$i H (k,i]) ps D'(k,r])P$ H (k). I then have the following 
result: 



o n 2 

cr LRS ~ fi67T 2 i ■ " 



-e 



o 12 o 



£ I k 2 dkP&(k,rio)P&(k,r)Q-£/k). (69) 



Now I argue that one can approximate the result by using £ ~ krjQ in the first power spectrum and 
using that P& (k) oc k~ 1 P(k) (which excludes the cancellation of (Mr/) -1 ) to give on small scales 
P& oc k~ 2 . Then it can be conveniently shown that: 

roo 

C„ NLRS ~ a (327r 2 ) e - 2K *£ / dkP# {k, % - l/k). (70) 

Jk H 

Here ao = 2A \&n / ?>Hq£Iq\ 2 . We once again use that k £(rjo — rf) and that r ~ (770 — rf) to then 
find that 

Cf^ s ~ £ 2 (32^ 2 ) f m %P^ {i y ri). (71) 

JO f T 

This can be readily compared to the Rees-Sciama calculation (where we ignore the damping): 

Nonlinear Rees - Sciama : C £ NLRS ~ £ 2 Cf s . (72) 

This is of course a gross approximation. However, I think that it does go some distance to demon- 
strate that if Cf s dominates the angular correlation functions near £ ~ 5000, one can expect the 
nonlinear Rees-Sciama effect (Eq. |5|) to dominate on scales much larger than this. 

1. If Cf s gives AT HS /T ~ 10~ 7 - 10~ 6 between £ ~ 100 - 300 (where the Rees-Sciama effect 



peaks |30|1), one can naively expect, from (72) to find that A2^ LRS /T ~ 10 5 — 10 4 between 
£ ~ 100 - 300. 

2. If I include damping and cancellation as k/k I then find : 

C< LKS ~£Cf s . (73) 

This then gives Aj^, lrs /T ~ 10~ 6 — 10~ 5 between £ ~ 100 — 300. The latter is probably more 
realistic. The keypoint is that it is does not seem to be a negligible effect. 



5 Conclusions 

Although my results on nonlinearity are by no means conclusive, I have evidence that the Newtonian 



threading suppresses nonperturbative nonlinearity on small scales, Eq. (44 - |50|). This along with 
the well know result that the Newtonian threading is inconsistent beyond linear order give an 
indication that such treatments are generically inadequate for the study of relativistic cosmology 
unless the a priori assumption that the universe is close to almost-FLRW on all observationally 



relevant scales is made. Although the self-consistency of this assumption has be shown [31, |3^] - 
such an approach is not generic. 

If this assumption is relaxed one discovers that there could be additional small-scale effects that 
have been excluded, by definition, from the canonical treatments of CMB temperature anisotropics 
Eq. (H) and Eq. @. One such small-scale effect, the Rees-Sciama effect, when included in the 
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frame work of the nonperturbative small-scale corrections to the radiation dynamics, leads one to 
the conclusion that there could be contributions of the same order of magnitude near the peak in 
the radiation transfer function as the anisotropies themselves due a coupling of the radiation via 
gravity to the nonlinear matter dynamics, Eq. (72). This conclusion is not consistent with the 
canonical treatment - from which it is excluded by construction. 

I have provided some evidence demonstrating that gravitational nonlinearity may become prob- 
lematic before the standard I = 5000 limit in the case of no feedback between the CMB temperature 
anisotropies - the scaling is different from the canonical treatment as an additional k scaling is in- 
troduced from the exact treatment. It can be expected, when the feedback is included, that the 
effect can become more problematic if there is no significant smoothing effect, Eq. (]73|). 

One weakness in the above approach is that the almost-FLRW anisotropies for high — I are 
found in the Newtonian frame, then boosted to the total frame. As ta ( is invariant under small 
frame boost the solutions are the same as those found in the Newtonian frame. The relationship 
between the shear and the acceleration potentials, the potentials and the matter perturbations and 
the peculiar velocities and the matter perturbations are all found in the total frame. Such tricks are 
only consistent for small relative velocities - at best this treatment is valid for mild nonlinearity. 
The other obvious weakness lies in the applicability of the string of approximations required to 
reduce the angular correlation function to that which takes on a form similar to the usual Rees- 
Sciama effect - which is the trick that allows the comparison between the Rees-Sciama effect and 
the nonlinear (in the Boltzmann equations) Rees-Sciama effect discussed here. 
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